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Abstract 

To  predict  random  fatigue  life,  the  popular  cumulative  damage 
criterion,  which  Is  based  on  constant  amplitude  sinusoidal  fatigue 
tests.  Is  known  to  be  inaccurate.  Therefore,  actual  random  fatigue 
experiments  are  proposed  to  establish  new  and  reliable  damage  criterion. 
In  this  report,  all  conceivable  probabilistic  factors  which  affect  random 
fatigue  life  are  explored  and  mathematically  analyzed.  Based  on  these 
factors  the  random  fatigue  experiments  will  be  planned  and  conducted. 
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INTRODUCTION 

There  are  two  major  approaches  for  Investigating  random  fatigue:  the 
cumulative  damage  approach  and  the  phenomeloglcal  approach.  The  main  concept 
of  cumulative  damage  criterion,  proposed  by  Miner  [1]  Is  based  on  the  results 
of  constant  amplitude  slnuslodal  fatigue  tests.  His  hypothesis  Is  charac- 
terized by  the  following  equation 

Z (n  /N  ) < c 
1 

In  which  n^  represents  the  actual  number  of  cycles  at  a given  stress 
level,  N^  the  fatigue  life  at  that  same  stress  level,  and  c Is  a constant. 
Miner  proposed  that  c should  have  a value  of  1.0.  It  Is  shown  [2]  that  the 
limited  experimental  results  which  have  been  used  to  evaluate  this  hypothesis 
show  a range  for  c from  0.3  to  about  3.0.  Since  sinusoidal  fatigue  results 
are  readily  available  this  criterion  became  very  popular  even  though  It  was 
repeatedly  proved  [3-6]  Inaccurate  to  predict  failure. 

The  phenomeloglcal  approach  Involves  actual  random  fatigue  experiments. 
Some  of  them  use  statistical  experiment  design  for  both  deterministic  and 
random  factors  and  most  of  them  use  statistical  and  variance  analyses  of 
results  [7-11].  The  present  report  explores  all  conceivable  random  factors 
which  affect  random  fatigue  life.  These  factors  are  then  mathematically 
analyzed.  The  mathematical  expressions  of  these  factors  facilitate  the 
design  of  random  excitation  signals  and  data  processing  of  random  responses. 
The  final  purpose  Is  to  establish  new  and  reliable  damage  criteria  for 
random  fatigue. 

ELASTIC  AND  PLASTIC  RANDOM  FATIGUE 

The  stress-strain  curves  of  tensile  tests  and  the  SN-curves  of  fatigue 
tests  Indicate  three  stress  levels  that  may  be  significant  to  random  fatigue. 
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They  are  fatigue  run-out  stress,  the  yield  point  stress,  and  ultimate  stress, 

represented  by  a^,  o^,  and  respectively.  The  corresponding  strain 

levels  are  denoted  by  e^,  e^,  and  e^.  The  yield  point  divides  the  elastic 

and  plastic  regions.  The  ultimate  strea,  not  rupture  stress.  Is  used  for 

both  brittle  and  ductile  materials.  In  the  case  of  brittle  materials  the 

ultimate  and  rupture  stresses  are  the  same.  Although  ductile  materials 

have  distinct  ultimate  and  rupture  stresses,  they  behave  like  brittle 

materials  under  dynamic  loadings.  Unless  all  the  peaks  of  strain  signals 

of  random  response  are  under  c^,  three  types  of  failure  problems  can  occur. 

The  first  two  types  are  random  fatigue  problems. 

If  all  the  peaks  of  random  strain  signals  are  below  c^,  and  there  are 

peaks  between  and  e^,  the  corresponding  fatigue  Is  elastic  random 

fatigue.  If  all  the  peaks  are  below  e , and  there  are  peaks  between  c 

u y 

and  e^,  the  corresponding  fatigue  Is  plastic  random  fatigue.  Once  the 
strain  Is  beyond  e^,  this  Is  a single  dangerous  state  of  first  excursion 
failure. 


PROBABILISTIC  FACTORS 

Based  on  the  discussion  of  elastic  and  plastic  random  fatigue, 
probabilistic  factors  conceived  to  affect  elastic  and  plastic  fatigue  are 
listed  as  follows. 

(1)  mean  strain 

(2)  variance  of  strain 

(3)  zero  strain  upcrossings 

(4)  and  level  upcrossings 

(5)  duration  of  excursion  between  zero  crossings 

(6)  duration  of  excursion  beyond  e,  and  e levels 

I y 
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(7)  peak  strain  probability  density  functions  - band  width 

(8)  average  peak  strain  amplitude  beyond  and  levels. 

MATHEMATICAL  ANALYSIS  OF  PROBABILISTIC  FACTORS 

For  each  random  signal  the  probabilistic  factors  will  be  mathematically 
analyzed  prior  to  the  random  fatigue  experiments.  The  statistical  expres- 
sions of  these  factors  other  than  mean  and  variance  are  derived  according 
to  Rice  [12].  The  derived  equations  are  general  and  applicable  to  random 
variables  and  random  processes  under  the  conditions  prescribed  In  the 
derivation;  and  they  will  be  applied  to  our  random  strain  signals. 

1.  Mean  Strain 

The  mean  or  expected  value  of  a continuous  variable  X of  a random 


process  X(t),  expressed  as  E[X],  Is  defined  as 

r ” 


E[X] 


xp^(x)  dx 


(1) 


where  the  probability  density  function.  It  is  assumed  that  the 

Integral  does  exist.  In  other  words  |x|p  (x)  may  be  Integrated  In  the 

A 

Interval  E[X]  Is  also  denoted  by  The  mean  value  represents 

the  steady  part  of  a random  signal  which  consists  of  steady  and  fluctuating 
components. 

2.  Variance  of  Strain 

The  expected  value  of  the  square  of  the  difference  between  random 
variable  X and  mean  y is  defined  as  variance  which  is  denoted  by 


a^.  Therefore 


aj  - E[(X-yj^)^] 


The  above  equation,  after  expansion,  can  be  written  as 


aj  - mh  - 


(2) 
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In  the  derivation  of  equation  (2)  the  following  expression 

r OO 

Pjj(x)  dx  « 1 

Is  used.  The  variance  Is  a measure  of  dispersion  and  Is  always  greater  than 

zero.  The  dispersion  can  be  measured  by  taking  the  positive  square  root  of 
2 

a^.  is  referred  to  as  the  standard  deviation,  or  root  mean  square  value 

(r .m. s . ) . 

3.  Zero  Strain  Upcrosslngs 

Strain  reversal  Is  considered  as  an  Important  factor  which  affects 
fatigue,  vniereas  In  sinusoidal  fatigue  It  Is  represented  by  the  number  of 
cycles.  In  random  fatigue  It  Is  represented  by  the  number  of  zero  strain 
crossings. 

To  Investigate  the  expected  number  of  crossings  of  a Gaussian  random 
process  x(t)  with  realizations  Xj^(t)  at  an  arbitrary  level  Xq  In  the 
time  Interval  At  - t2-tj^,  it  Is  necessary  to  construct  a counting  functional 
n.  The  unit  step  function  u(t)  Is  used  for  this  purpose. 

Let  a realization  of  a new  random  process  Y(t)  be  defined  as 

Yj^(t)  - u{xj^(t)-XQ} 

Differentiating  with  respect  to  time  gives 

Yj^(t)  - ij^(t)  d{xj^(t)-XQ} 

which  vanishes  everywhere  except  when  Xj^(t)  - x^,  at  which  point  there 
exists  a spike  of  unit  area  directed  upward  or  downward  depending  on  whether 
ij^(t)  Is  positive  or  negative.  The  counting  functional,  i.e.,  the  nuiid>er 
of  crossings  per  unit  time  at  t,  for  this  typical  realization  of  the  random 
process  x(t),  given  by  Middleton  Il3]  as 

n(xQ,t)  - |xj(t)|  6{x^(t)-XQ} 
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Includes  both  upward  and  downward  crossings.  The  total  number  of  crossings 


In  the  Interval  (t^,t2)  Is  therefore 


N(Xo.ti,t2) 


r t. 


J •'i 


|Xj^(t)|  6{xj^(t)-XQ}  dt 


The  expected  number  of  crossings  In  this  typical  realization  x^(t)  of  the 
Gaussian  random  process  x(t)  Is  given  by 

r 

E[N(xQ,t^,t2)]  " J =Ui(t)|  6{x(t)-XQ}]dt 


r ^2  r “ r 

J J —00  J — oo 

(tj-tl)  j 


x|  6(x-Xq)p^(x,x)  dx  dx  dt 


where  p^(x,x)  Is  the  joint  probability  density  function.  Since  each 
upcrosslng  Is  followed  by  a downcrossing  the  expected  nuiid>er  of  upcrosslngs 
per  unit  time  Is  given  by 

^ E[N(xQ,tj^,t2)] 


2 >-1 


1 

2 


From  the  definition  of  probability  density  function,  we  have 

1 


Pxi(***)  “ — 

2n.^  m-(l-pZ) 


0 2 


exp  i-^  (^-2ei-.si), 
2(l-p^  "o  ,5^  *2 


In  which 


_ ELx(t)i(t)] 

/inr 

^0*2 


i 


where  o)  is  the  angular  frequency  and  (hi)  Is  the  one  sided  spectral 


density  function. 

For  ergordlc  process  the  expectations  are  equal  to  the  corresponding 
temporal  averages.  Thus  ] 

E[x(t)x(t-K)]  « Rjqj(t)  ■ <x(t)x(t+T)>  i 

where  R__(t)  Is  the  autocorrelation  function  and  < > Indicates  the  1 ' 

temporal  average.  We  also  have  ; ! 


and  the  equation  for  E[N^.(Xq)]  becomes 


0 the  expected  number  of  zero  upcrosslngs  per  unit  time  Is 


obtained  as 


4.  e,  and  e Level  Upcrosslngs 


The  expected  number  of  upcrosslngs  at  the  level  x.,  as  derived  above 


Therefore  the  equations  for  the  expected  number  of  upcrosslngs  of  a random 
strain  signal  at  the  levels  and  per  unit  time  are  obtained, 

respectively,  as  follows. 


5.  Duration  of  Excursion  Between  Zero  Cro8s<nffa 


The  duration  of  the  excursion  Is  the  time  Interval  during  which  the 
oscillation  exceeds  a specific  level.  A particular  case  Is  the  duration 
of  excursion  between  zero  crossings. 

Let  an  arbitrary  level  x^  be  exceeded  by  the  realization  Xj(t)  of 
the  stationary  random  process  x(t)  at  an  upcrosslng  when  the  oscillation 
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has  positive  slope  x(t^)  > 0 while  at  a downcrossing  the  oscillation  has 
negative  slope  xCt^)  < 0.  Define  a realization  of  a new  random  process  as 

^1^^1’^2^  “ u{xj^(tj^)-XQ}u{ij^(t^)-0}u{xj^(t2)-XQ}u{-*j^(t2)-0} 

In  which  u Is  a step  function.  Differentiating  with  respect 

to  t^  and  then  with  respect  to  t2  gives  four  terms.  One  of  these  terms 
Is 

Xi(ti)6{Xi(ti)-Xo}u{Xi(ti)-0}Xi(t2)6{xi(t2)-Xo}u{-Xi(t2)-0} 

which  represents  a spike  of  unit  area  when  Xj^(tj^)  = “ *0’  *l^*’l^^  ® 

and  x^(t2)  < 0.  It  can  be  used  as  the  counting  functional  representing  the 
number  of  crossings  per  unit  time  whenever  ■ Xj^(t2)  “ x^  and 

provided  that  ^ ® ^ this  counting  functional 

be  denoted  as  n(xQ,XQ,Xj^,X2,t) . Then 


n(Xo,Xo,Xi,X2,t)  = Xj^(tj^)Xj^(t2)6{xj^(tj^)-XQ}u{xj^(tj^)-0} 
6{xi(t2)-Xo}u{-3ti(t2)-0} 

For  a zero  mean  valued  Gaussian  random  process,  let  t^^  * t and  t2  “ t + T, 
the  total  number  of  crossings  for  duration  of  excursion  t at  level 
during  time  interval  ^2  ~ ‘"I 


N(Xo»*o’*’*t’^1’*^2^ 


n(xQ,XQ,x,x^,t)  dt 


where  x » x(t),  x^  * i(t+T).  The  expected  number  of  the  duration  of 
excursion  at  level  x^  for  the  time  Interval  t2  - tj^  or  x Is 

E[N(xQ,XQ,x,x^,tj^,t2)] 


r t. 


E[x^(tj)Xj^(t2)6{x^(tj^)-XQ)u{xj^(tj^)-0}6{xj^(t2)-XQ}u{-x^(t2)-0}l  dt 


10. 


-(tj-tj) 


0 r <*> 


* Px^*0’*0’***T^  *** 


J -oo  J 0 


The  expected  number  of  excursion  Interval  t per  unit  time  Is 

0 /■  » 


E[N(T)]  - - 


In  the  above  equation  the  four  variable  density  function  Is  given  by 

-axJ/|A^I 


1 1/2  e " ” exp[-  [bx^Cx-x^)  + dxx, 

. e /•2  . 2. , 

+ 2 ^ 

where  Is  a four  variable  variance-covariance  matrix,  expressed  as 


“o 

mo(T) 

0 

nij^(T) 

mo(T) 

“o 

0 

0 

-nij^(T) 

“2 

02 

0 

012  (t) 

“2 

and 


a - {m2-m2(T)}t{mQ-mQ(T)}{m2-hn2(T)}  - m^(T)] 
b - “i(T)[{mQ-mQ(T)}{m2-hn2(T)}  - mJ(T)] 
d - -m2(T){mJ-mJ(T)}  + mQ(T)mJ(T) 

e - ■ “o“l^‘'^^ 

■o(T)  - Etx(t)x(t+T)]  - 
■j(T)  - Elx(t)x(t+T)]  - R^(t) 

■2(t)  - E[x(t)x(t+T)]  - R..(t) 


The  determinant  of  is  expanded  and  simplified  as 

l\l  “ [{niQ+mQ(T)}{m2-m2(T)}  - m^(T)  ] [ {mQ-mpCx) }{m2+m2(T) } - m^(T)] 

2 ^2 
e -d 


m2-m2(x) 


The  probability  density  function  for  the  excursion  interval  X between 

= 0 ii 
E[N(x)], 


zero  crossings  for  which  x^  = 0 is  defined  as 


v® 

E[N  (0)] 


where 


1 - ii 


E[N(X)] 


X '0  .21 

o 1 


noo 

**x®^p  li;;!  ^‘^“x  2 


Rice  has  evaluated  the  above  Integral  for  E[N(x)]  _ and  gives 

Xq-O 

Po^”"^  " [1  + h cot-\-h)] 


where 


h - d(e^-d^)  ^ 


Therefore 


E[x] 


Xq-O 


r T, 


xpQ(x)dx 


max 


xpQ(x)dx  (5) 


6.  Duration  of  Excursion  Beyond  and  Levels 

The  probability  density  function  for  excursion  Interval  x between 
an  upcrosslng  and  a downcrossing  at  level  x^  is  defined  as 

p(x)  . — 

E[N^.(Xq)] 


12. 


where 


EFN 


E[N(t)] 


-(Xo/2mo) 


as  obtained  previously.  Instead  of  a tedious  numerical  Integration  of  the 
above  equation  for  E[N(t)],  Tikhonov  [14]  derives  an  approximate  expression 
for  the  probability  density  function  for  the  duration  of  excursion  at  an 
arbitrary  level  Xq  as 

P(t)  - pQ(T)exp(-a^/4) 


where 


Then 


E[N(t)] 


J 


r T 


Tp(T)dT 


max 


Tp(T)dT 


(6) 


7.  Peak  Strain  Probability  Density  Function  - Band  Width 

A maximum  of  a typical  realization  Xj^(t)  of  the  Gaussian  random  process 
x(t)  occurs  provided  that  the  realization  x^(t)  of  the  random  process  x(t) 
Is  zero  and  the  realization  x^(t)  of  the  random  process  x(t)  Is  negative, 
l.e.  Xj^(t)  - 0,  Xj^(t)  < 0. 

Define  the  counting  functional  as 

n(xQ,0,t)  - -Xj^(t)  6{ij^(t)  - 0}u{xj^(t)  - x^} 

tdilch  represents  a spike  of  unit  area  whenever  x^(t)  > x^,  i^(t)  ■ 0,  and 
x^(t)  < 0.  The  total  number  of  peaks  In  the  time  Interval  At  ■ t2  - tj^  Is 


13. 


N(Xq»0, tj^, t2)  “ ” 


2 
J 


Xj^(t)  6{xj^(t)  - 0}u{x^(t)  - Xq} 
t- 


dt 


Then  the  expected  number  of  peaks  in  this  typical  realization  x^Ct)  of  the 
Gaussian  process  x(t)  is  given  by 

f 

E[N(xQ,0,tj^,t2)]  - - E[xj^(t)  6{xj^(t)  - 0}u{xj^(t)  - x^}]  dt 

^ *^1 

The  expected  ntiiid>er  of  peaks  greater  than  x-  per  unit  time  for  a zero 


mean  valued  Gaussian  random  process  x(t)  is 

00  f ^ r 00 


E[N(Xq,0)] 


dx 
-00  J 


dx 


X(S{^0}u{x-Xq}p  ...(x.x.x)  dx 


dx 


) *0 


...(x.OiX)  dx 


XXX 


In  the  above  equation,  the  probability  density  function  is  derived  as 

1 ^ 1 , 2 . » 2 .. 


p ...(x,0,x) 

XXX 


from  which 


15. 


PY(y) 


^Pjjjyj(y.O,x)  dx 


Substituting  the  probability  density  function  and  Integrating  gives 


the  following  probability  density  function 

-yj/rm^f2 

e 


PY(y) 


-y2/2moe"  ^ 


-y^/2ni„ 


y e 


[0.5  + 


0 


where 


1 - 


”o“4 


r y 


erf(y)  - -==■ 

H 211 

From  Integration  table 


J 


-z^n 

e dz 


[ ° 1 

1 

f 0 

erf  i 

- 1 

1“ 

1 

[0.5j 

The  quantity  e which  can  be  used  as  a factor  to  classify  the  probability 
density  functions  and  evaluated  by 


in« 


1 - 


V4 


(7) 


Is  referred  to  as  the  band  width. 

8.  Average  Peak  Amplitude  Beyond  e.  and  e Levels 
-1 Z 

Let  N be  the  total  number  of  signal  peaks  In  a given  time  Interval 
and  Uq  be  the  number  of  peaks  with  magnitude  greater  than  a given  level 
x^  In  the  same  time  Interval.  Then  the  number  n Is  defined  as 


N 
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The  probability  that  the  peaks  exceed  an  amplitude  x.  Is  given  by 


*0^  “ n “ ^ 

J 

As  the  peak  distribution  of  a random  process  Is  described  by  the  Rayleigh 

probability  density  function,  so  that 

r » 2 . 

— * — e dx 

n m. 

X*  0 

J o 


from  which  n Is  obtained  after  Integration  as 

XQ/2mQ 


The  average  value  of  the  — th  highest  peaks,  denoted  as  x.  Is 

n 1 


expressed  as 


n ^ 0 

^ 00 

J *0  "" 


xp^(x)dx 


-x^/2m^ 


After  Integration 


■ ”/^o  (An  n)^  + /fiT  [0.5-erf(2An  n)^]} 


where 


erf  X 


r 

, 1 

ifTv 

J 


* V/2 


Finally 


x^/2bi  ^ 

*l“*  ° ° +//^  {0.5-erf (xJ/bq)^}]  (8) 

n 

for  average  peak  amplitude  beyond  e,  and  e levels  provided  x.  is 
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equal  to  and  e^,  respectively. 

CONCLUSION 

After  defining  elastic  and  plastic  random  fatigue  In  this  report,  all 
probabilistic  factors  which  are  conceived  to  affect  random  fatigue  are 
analytically  studied.  In  addition  to  the  mathematical  definitions  of  mean 
and  variance  all  other  factors  are  mathematically  analyzed  and  expressed 
In  terms  of  moments  of  spectral  density  functions.  Based  on  some  or  all  of 
these  probabilistic  factors  random  fatigue  experiments  will  be  planned  and 
conducted . 
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